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ABSTRACT 

A detailed gravimetric geoid of North America, the North Atlantic, Eurasia, 
and Australia computed from a combination of satellite -derived and surface l°x 1° 
gravity data, is presented herein. Using a consistent set of parameters, this geoid 
is referenced to an absolute datum. The precision of this detailed geoid is ±2 meters 
in the continents but may be in the range of 5 to 7 meters in those areas where data 
was sparse. Comparisons of the detailed gravimetric geoid with results of Rice for 
the United States, Bomford and Fischer in Eurasia, and Mather in Australia are pre- 
sented. Comparisons are also presented with geoid heights from satellite solutions 
for geocentric station coordinates in North America, the Caribbean, Europe, and 
Australia. 
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1. INTRODUCTION 

This paper presents a detailed gravimetric geoid of North America, the 
North Atlantic, Eurasia and Australia based upon a combination of satellite de- 
rived and surface 1° x 1° gravity data. Early gravimetric geoid computations 
were based upon surface gravity data; Hirvonen (1934) and Tanni (1948, 1949). 
The most ambitious of the pre-satellite gravimetric geoids was the Columbus 
geoid (Heiskanen, 1957). All of these pre-satellite geoids suffered from a lack 
of worldwide gravity coverage. With the advent of satellites it has been possible 
to derive the long wavelength components of the gravity field on a worldwide 
basis with considerable accuracy. The satellite gravity data can be combined 

with surface gravity data, in areas where surface gravity is available, to pro- 
vide accurate estimates of the details of the geoidal undulations. This ability 
to combine the two-data types to obtain detailed geoid undulations in local areas 
through combination of surface and satellite gravity was recognized early in 
satellite geodesy (Khan and Strange, 1966) and awaited only the gathering of 
sufficiently accurate data. The method of Khan and Strange (1966) has been 
applied essentially unchanged to derive the results presented in this paper. 

The geoid is becoming increasingly important for the support of research in 
geodesy and geophysics. The Skylab and GEOS-C spacecraft will carry radar 
altimeters for the purpose of measuring the geoid undulations in oceanic areas. 
An independently derived geoid map will provide a valuable complement to these 
experiments. By studying the gravimetric geoid, optimum experiment locations 
can be established. Also, the gravimetric geoid can be used to calibrate and 
check the accuracy of the altimeter. An accurate geoid map is also valuable for 
satellite and inertial navigation systems which are being used for offshore min- 
eral exploration. A number of experimenters have derived values for tracking 
station coordinates from satellite observations. For the cases where the height 
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of the station above mean sea level is known accurately, the geoid map can be 
used to check the accuracy of the derived height above the ellipsoid in past solu- 
tions. Furthermore, the detailed geoid can be used as a constraint for future 
solutions as was recently done by Mueller and Whiting, 1972 who incorporated 
an earlier GSFC detailed gravimetric geoid map (Vincent et al. , 1971) into their 
global geometric solution. Over the continents the gravimetric geoid can be 
compared with astrogeodetic geoids. These comparisons provide not only an 
indication of the relative accuracy of the geoids but they also provide information 
on tilts of the major geodetic datums with respect to a geocentric reference 
system. 

In previous publications (Strange, et al. , 1971; Vincent, et al. , 1971) de- 
tailed geoid height maps were presented covering a substantial part of the 
northern hemisphere. As was stated in the previous papers the geoid heights 
given were only accurate to within an additive constant. This constant difference 
between the geoid heights presented and true geoid heights included a zero order 
undulation error (Rapp, 1967) plus other small systematic errors. Since the 
presentation of these previous geoids, a study has been made to determine best 
estimates of the various parameters which could cause the systematic errors 
(Strange and Richardson, 1972). Using a consistent set of parameters based 
on the results of this study it is now possible to remove the systematic error 
to a large extent. The objective of this paper is to present a set of detailed 
gravimetric geoid height maps which are based on a consistent set of parameters 
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that are very nearly true values for the earth. Systematic errors in the geo- 
centric radius vectors to the geoid computed using the geoid heights of these 
maps and a reference ellipsoid of a^ = 6378. 142kms and f = 1/298.258 are al- 
most certainly less than 5 meters, i. e. an accuracy for the radius vector of 
better than lppm. 

In addition to being referenced to an absolute datum, the maps presented 
here also differ from those previously published in that additional observational 
data has been used and the area covered in previous reports has been increased 
to include Australia and Canada. 

The detailed gravimetric geoids presented here have a precision of ±2m 
rms. This precision was established by comparing the detailed gravimetric 
geoid with Rice's astrogeodetic geoid for the United States, Bomford’s astro- 
geodetic geoid for Europe and Mather's astrogeodetic geoid for Australia. 

Comparisons have also been made between the detailed gravimetric geoid 
and satellite derived tracking station positions of Goddard Space Flight Center 
(GSFC) and Smithsonian Astrophysical Observatory (SAO). 

2. METHOD OF COMPUTATION 

The geoidal undulation at any point P on the earth can be computed using 
the well known Stokes' formula: 

2TT tt/2 

N «- X >= 4 n G ' / /. A g T (T^ *0 S(y) cos cp'dq/dX' (1) 

x '=° 
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where: 


cp, X. - The geocentric latitude and longitude, respectively, of the computation 
point. 

cp , A = The geocentric latitude and longitude, respectively, of the variable 
integration point. 

N(cp, X) = Geoid undulation at cp, X . 

R = Mean radius of the earth. 

G = Mean value of gravity over the earth 
^g T ( c P*^ ) = Free air gravity anomaly at the variable point 9 ', X'. 

S(Y) = 'iT n(Y/2 ) “ 6 sin ^ /2 > + 1 + 5 cos Y 

- 3 cos Y In (sin OF/2) + sin 2 (Y/2)) 

where 

Y = cos^Csin cpsin cp'+ cos cpcos cp'cos(A - /)] 

In order to combine surface and satellite gravity data for geoid computation, 
the earth is divided into two areas, a local area (A-j) surrrounding the point P, 
and the remainder of the earth (A^) . Also the anomalous gravity in each area 
is partitioned into two parts represented by the symbols Ag and Ag„. The Ag 
values are defined as that part of the anomalous gravity field which can be repre- 
sented by the coefficients in a satellite derived spherical harmonic expansion of the 

gravitation potential. The 1969 SAO Standard Earth (Gaposchkin and Lambeck, 1970) 
was used in all computations described in this paper. The Ag values are defined as 
the remainder of the anomalous gravity field. Using this division of the earth's surface 
into two areas and of the anomalous gravity into two components one can write equa- 
tion ( 1 ) in the form: 
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( 2 ) 


N(cp, X) = Nj + N 2 + N 3 


where 


TT 


N = -5- 
1 4nG 


/ / 

0 TT 

~~2 


A * / 


Ag (9« X ) S(Y) cos cp dcp dX 
s 


N. 


R 


2 4ttG 


// 


Ag (cp \X') S(Y) cos cp'dcp'dX' 

Li 


(3) 


N = JL 

3 4ttG 


u 


Ag. (cp*, X*) S(Y) cos cp' dcp ” dX 
z 


The following paragraphs discuss how each of the three components presented 
in equation (3) is handled in the computations. 

Given a set of satellite derived coefficients in the spherical harmonic expansion 
of the gravitational potential, a number of methods exist for the computation of the 
Ni component of the geoid undulation. 

The computation of was not carried out in the present case by using the 
integration indicated in equation (3). Rather the procedure described by Bacon, 

et al. , (1970), was used. Briefly this procedure consists of fixing a value of the 
potential, W q , and computing the component N as 

Ni = r - r E (4) 
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where: 


r = is the radial distance to the equipotential surface defined by 

W q and the potential coefficients of the SAO 1969 Standard Earth. 

r„ = is the radial distance to a selected reference ellipsoid defined by 

hi 

a semimajor axis (a g ) and flattening (f). 


The radial distance, R^, to the equipotential surface W q at a particular latitude 
and longitude 0 , is determined by using the equation 




where 

kM = the product of the gravitational constant and the mass of 
the earth 

a e = semimajor axis of the reference ellipsoid 
r = geocentric radius 
cj = earth's angular velocity 

C nm and S nm = fully normalized spherical harmonic coefficients of the 
gravitational potential 

P nm (sin0) = the Associated Legendre Polynomial. 

The only unknown in this equation is r. Values of r-^ = R + e , r 2 = R - e, and 
r 3 = ( r l + r 2 )/^ are chosen for substitution into equation (5) for evaluation of the 
functions 

0 1 (r x , 0,\),0 2 ( r 2 » (p. A.), and 0 g (r g , 

The r- for which U. - W I is a maximum is identified and eliminated from 
1 1 o 

consideration. The two remaining values of r^ are labeled r-^ and r 2 and are 
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used for calculation of r 3 = (r± + r 2 )/2. The potential functions are evaluated with 

\ 

these arguments and the worse-value elimination process is repeated. The process 

-12 

continues until an r is chosen such that \\p (r, 0, k) - W | £ 10 • Using this 

value of r and the value of r^ computed using the input values of a e and f of the 
reference ellipsoid, a geoid undulation component N-j_ is computed. 


For the computations described in this paper, the area for a point at 
which the geoid was being computed was defined to consist of a twenty degree by 
twenty degree area centered on the computation point. The computational formula 
used was: 

400 _ 

N 2 = iFFG" 23 % (9.', Vj) S(Tj) cos cpj Acp'AX' (6) 


where 


Ag_(9 * , is the mean value of Ag within the 1° x 1° square 
2 J 1 2 

th 

S(Y .) is the value of Stokes' function at the center of the j 1° x 1° square. 

Acf/- hk' = 1°. ' 

The value of Ag used for each 1° x 1° square was computed using the formula 


ig 2 = 4g e - \ 


( 6 ) 


data. 


The Ag^ values are mean 1° x 1° free-air anomalies provided by surface gravity 
The primary source material for this gravity data is as follows: 


1) North America: Strange and Woollard (1964), Woollard (1968), Nagy (1970) 
and ACIC (1971) 


2) North Atlantic: Bowin (1971), Talwani (1971), Strang Van Hees (1970), 
and ACIC (1971) 
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3) Eurasia: Tengstrom (1965), Arnold (1964). Bowin (1971), and ACIC (1971) 

4) Australia: Mather (1969). 

Values of Ag e for each 1° x 1° square were computed by carrying out the computation 

Ag = Ag + Y - P.C. - Y 
s e S IF IF N 

where 

Ag T „ = mean value of free air anomaly referred to the international 
IF 

gravity formula 

Yj F = international gravity formula 

P.C. = Potsdam correction with a value of -13.7 mgals 

Y^ = 978.0322 (1+. 0053025 sin^ cp- .00000585 sin^ 2cp) gals 

In carrying out the computations Y^, and Y^. were evaluated at the center of each 

,0 ,o 

lxl square. 

The Ag s values are that part of the mean 1° x 1° free-air anomalies represented 

by the satellite harmonic coefficients used in computing N,. The Ag values are 

l s 

obtained by evaluating the following equation at the center of each 1° x 1° square. 

— k n _ _ 

Ag = y y y (n-1) [ C cos ml'+ S sin mX' ] P (sind,') (7) 

s e nm nm nm ^ ' ' ; 

n=2 m=0 

where 

Y = Equatorial gravity in milligals 
© 

k = Upper limit on degree and order of the geopotential model 

n = Degree index of harmonic coefficients 
m = Order index of harmonic coefficients 
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In equation (7) , the C and C terms do not represent the complete coefficients but 

mU tcU 

rather the difference between the complete coefficients and the coefficients compatible with 
the ellipsoid used in computing N^. The difference values used were AC^q = .03577 x 10 ^ 
andAC 4Q = -.232 x Kf 3 * * 6 (fully normalized). In order for the above described proce- 
dure to produce correct results, the quantities Ag , Ag , and the a and f which define 

the ellipsoid used to compute must all be compatible. Compatibility implies that 
the values of and used to compute the values of theoretical gravity needed to 
obtain Kg and Ag are the same as the values of C and C implied by the reference 

G S ZU tU 

ellipsoid. Correct results in the absolute sense are also dependent upon the value of 
W q chosen to represent the true value of the potential of the geoid. The effects of 
not making Ag , Ag , a, and f compatible are twofold. First, all the computed geoid 

G S 

heights may by in error by a constant; in addition, there will be a systematic error as 
a function of latitude. The effect of selecting an incorrect value of would be to 
introduce a constant error in all geoid heights. 

In the calculations described here the term N in equation (2) is set equal to zero. 

O 

this is equivalent to assuming that the satellite derived approximation to the gravity 
field is adequate for the area at a distance of greater than ten degrees from the 
computation point. 

3. GEOIDAL SCALE 

With the procedure described in Section 2, one can compute values of geoid 
height, N, which are accurate to within an additive constant provided that one; 1) uses 
the correct value of rotation rate, w, 2) assumes that the reference gravity formula 
used in computing the surface gravity anomalies is compatible with the flatten- 
ing chosen for the reference ellipsoid selected and 3) makes the A C and AC used 

Zv rrl) 

in deriving the satellite component of geoid undulation compatible with the selected 

reference ellipsoid flattening. With these conditions satisfied, one need only use 

reasonable approximations to such parameters as W q , kM, Potsdam correction, a^, 

and V to obtain values of N accurate to within an additive constant, 
e 
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It is desirable, however, to produce geoid heights which do not have a constant 
uncertainty and that the parameters of the best fitting reference ellipsoid be known. 

The best fitting reference ellipsoid is defined here as that ellipsoid having semimajor 
axis, a g , and flattening, f, such that its surface is an equipotential having the same 
value of the potential as the geoid, when the potential is computed using the kM, and 
J 2 of the actual earth. As a result of recent work (Strange and Richardson, 1972) it 
now appears possible to present absolute geoid heights and reference ellipsoid param- 
eters with reasonable assurance that any systematic errors in geocentric radius vector 
to the geoid are less than 5 meters and probably of the order of 3 meters. 

In using the theory described in Section 2. systematic errors can occur due to 

errors in the following parameters: kM, u>, J , a , f, Y , and the Potsdam cor- 

2 e e 

rection. Strange and Richardson (1972) adopted values of 

9fl q o 

kM a+ j= 3.986012 x 10 /cm /sec 

90 9 9 

kM = 3.986009 x 10 cm /see 
e 

» = .72921151467 x 10~ 4 rad/sec 

J 2 = 1082.6392 x 10~ 6 

1/f = 298.258 

Potsdam correction = -13.7 mgals 

taken from previous work together with available surface gravity values to compute 

Y g = 978. 0332 gals 

a = 6378. 1388 kms 

e 

W q = 6263691.0 kgal m 
where 

M 

e 
M 

a+e 


mass of the earth excluding the atmosphere 
mass of the earth including the atmosphere 



It was also found by comparison of detailed gravimetric geoid heights with 

dynamic station positions, that the values (after correcting the dynamic station positions 

to be compatible with the chosen values of kM ) of a ranged from 6378. 141 kms to 

a+6 e 

6378. 144 kms. If these values were correct, commensurate changes in Y and W 

e o 

would be required to maintain a consistent set of parameters. Since the value of 

Y g - 978.0332 gals obtained from surface gravity analysis could easily be in error 

by±l. 0 mgals, it was decided to adopt as the set of parameters for the geoid undulation 

the k M , kM , id, J and l/f given above together with 
a+e g 2 

Y = 978.0322 gals 

e 

a = 6378.142 kms 

e 

W = 6263687.5 kgal m 
o 

Geoid undulations were then computed using this coherent set of parameters. 

As a means of evaluating the scale of the geoid, detailed geoid heights and 

reference ellipsoid parameters were used together with mean sea level heights taken 

from the NASA Directory of Observation Station Locations (NASA, 1971) to compute 

geocentric radius vectors for a substantial number of satellite tracking stations . These 

geocentric radius vectors were then compared with the geocentric radius vectors 

derived from satellite observations by a number of investigators in order to determine 

the extent of systematic differences (tables 1 through 3). Since the scale of the dynamic 

orbit analysis is set by the value of kM used, the dynamically derived radius vec- 

a+e 20 3 2 

tors were modified to be compatible with kM = 3. 986012 x 10 cm /sec before 

3. ' G 

making any comparisons. The dynamic radius vectors and those obtained using the 
gravimetrically derived parameters have a systematic mean difference of 3 meters or 
less. This level of agreement must be considered excellent taking into account the 
potential uncertainties in the various data used in deriving the computational param- 
eters. Of the various potential sources of the differences, the most probable 
causes are: 
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1) Errors in detailed gravimetric geoid height at tracking stations 
due to the use of simple free air anomalies rather than 
terrain corrected free air anomalies. 

2) Errors in values of Y . W n and a_ 

GO c 

3) Errors in dynamic station coordinates including small deviations 
from the center of mass of the earth of the origin of the coordinate 
system to which dynamic station positions are referenced. 

4) Errors in mean sea level elevations for some tracking stations. 

Theoretically terrain corrected free-air anomalies rather than simple free air 
anomalies provide more accurate estimates of geoid height. The effect of using simple 
free air anomalies is to produce geoid heights which are systematically too negative in 
the vicinity of land areas with rugged relief. Dimitrijevich (1972) has shown that the 
value of the difference in the United States ranges from in excess of +3.5 meters in 
the rugged mountains of the western United States to about +0.2 meters along the east 
of the United States. Since most tracking stations used in the comparisons are on 
large land masses and several are in areas of rugged relief, one to two meters of the 
systematic difference can be assumed to arise from this source. It should be noted 
that differences due to this source are not the result of errors in basic parameters 
but the use of a slightly incorrect form of surface gravity anomalies in the computations. 

For reasons which are not entirely clear it has been found that dynamically 
derived station positions of different investigators can have small systematic differ- 
ences in the X, Y, Z values on the order of 5 meters in magnitude. It should be noted 
here that the GSFC long arc solution station coordinates have been modified to include 
a 10 meter correction to the Z values to account for a shift of the coordinate system 
from the center of mass. Because the tracking stations used in deriving the systemtic 
effects given in tables 1 through 3 are not uniformly distributed over the earth's sur- 
face one could anticipate the possibility of contributions of a few meters due to the 
failure of the origin of the coordinate system to coincide with the center of mass of the 
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earth. Of course the chosen value of a = 6378. 142 kms could still be in error by 

e 

1 to 3 meters. 


Thel causes \of the small differences noted in Tables 1 through 3 for the grav- 
imetric geoid heights and those obtained from the dynamic station positions of most 
investigators are no doubt the result of some combination of the above four error sources 
together with a large number of other small error contributors. However it should be 
noted that the above comparisons give no information concerning possible errors caused 
by the adopted values of k M^ +e and k M . It seems unlikely these would exceed two to 
three meters. 


4. PRESENTATION AND EVALUATION OF RESULTS 

A detailed gravimetric geoid, a detailed gravimetric geoid - satellite geoid, and 
a satellite geoid for North America, the North Atlantic, Eurasia, and Australia were 
computed. These are presented as Figures 8 through 16. 

To evaluate the precision of the detailed geoid for the areas computed, a number 
of comparisons were made. The first comparison was made with the astrogeodetic 
geoid data of Rice (1970) for the United States. Before any comparisons could be made, 
Rice's data were transformed from the North American Datum (NAD) to the geocentric 
coordinate system using several transformation sets. Table 4 presents the differences 
between Rice's astrogeodetic geoid and the gravimetric geoid using five different sets of 
translation elements after removing the mean differences. In all cases, the rms differ- 
ences are on the order of 2 meters or less . 

In Europe a comparison was made with Bomford's (1971) astrogeodetic geoid 
map. Bomford's astrogeodetic geoid values were first transformed from the European 
datum to the geocentric system using GSFC Long Arc Solution transformation sets of 
AX = 89, AY = 120, and AZ = 118 meters. It should be noted that these are mean 
translational values that do not incorporate the tilt of the European datum with respect 
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to the geocentric system. The comparisons were made along profile lines at latitudes 
44°, 48° N and longitude 9° E. The comparisons of the transformed astrogeodetic 
geoid and the gravimetric geoid along the latitude profiles show an east-west rotation 
in the European datum which is not shown in the longitude profile (see Figures 1 through 
3). This rotation was on the average equal to about 1.7 arc seconds. Allowing 

for this rotation in the datum, the relative agreement between Bomford’s 

astrogeodetic geoid and the gravimetric geoid is within the ±2 meter range. 

Another comparison in Eurasia was made with Fischer’s astrogeodetic geoid map. 
Table 5 shows a comparison of the detailed geoid and the results of Fischer (1968) for a 
traverse across Eurasia at 52° N latitude. Both sets of geoid heights are referred to 
an ellipsoid of a^ = 6378. 142 km, Fischer's values being taken directly from her 
paper. The geoid heights agree very well if one assumes a systematic error of 
13 meters for Fischer's values west of 50° E longitude and 21 meters for values 

east of 50° E longitude. The cause for this disagreement in Europe has not been 
explored but could easily arise from a 10 to 15 meters error in the values used by 
Fischer to transform to the geocentric system as well as by the fact that she used only 
translation parameters to transform to the geocentric system while our investiga- 
tions show that a rotation is also required. The change of 8 meters in the systematic 
difference at 50° E longitude could either arise from errors in connection of the 
Pulkova datum to the European datum or from systematic errors in the geophysically 
predicted free air anomalies used for gravimetric geoid computations in this area. 

In Australia comparisons were made with Mather's Astrogeodetic geoid (1971). 
Again, before any comparisons were made, Mather’s astrogeodetic data were transformed 
to the geocentric system using Mather’s transformation values of AX = -125, AY = -30, 

A Z = 111 meters. The comparisons were made along profile lines at latitudes 22°, 

26 , and 30 S (Figures 4 through 6). The comparisons along these latitude profiles 
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show an east-west rotation in the Australia datum approximately equal to 0. 5 
arc seconds. This rotation has also been noted by Mather (1971). Allowing for 
this rotation in the datum the relative agreement between Mather's astrogeodetic 
geoid and the gravimetric geoid is within ±2 meters. 

5. COMPARISON OF SATELLITE GEOID WITH DETAILED GEOID 

It is of interest to compare the satellite derived and detailed geoid maps. 
This comparison has been performed in two ways; by comparing latitude profiles 
and also by plotting contour maps of the differences. 

Figure 7 presents a profile drawn across approximately one-half the globe 
at latitude 35°N comparing the detailed gravimetric geoid and the SAO '69 satel- 
lite geoid. The SAO '69 satellite geoid was referenced to a e = 6378. 142km and 
a kM value of 3.986012 x 10 20 em 3 /sec 2 . Several conclusions result from the 
study of this profile: 

1) The existence of a steep gradient in the North Atlantic . 

2) There is no indication on the profile of any major tilt in North America. 

3) There is an east-west tilt in the geoid of Europe. 

4) In general the satellite geoid differs at most by approximately 10 m from 
the gravimetric geoid. 

Figures 9, 12, and 15 present contour maps of the differences between the 
detailed geoid and the satellite geoid for Australia, North America, the North 
Atlantic and Eurasia. In Australia the most prominent differences occur in 
areas east of 135° longitude where they range between 10 and 20 meters. These 
large differences are attributed to the dominance of the eastern mountain ranges 
which adjoin relatively flat plains on the west and a shallow continental slope 
on the east. Table 6a lists coordinates for twelve locations in North America 
and the North Atlantic where the difference is in excess of 8 meters. The 
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difference of -12 meters over the Puerto Rican Trench was not unexpected since 
the gravity gradient is large over a small region. Other areas, for example; 
42°N, 72°W and 49°N, 61°W where the differences are 14 meters and 10 meters 
respectively, may indicate broad shallow features to which satellites are not 
sensitive. Table 6b contains thirteen locations in Eurasia where the geoid height 
difference is larger than 8 meters. The relatively poorer agreement in Eurasia 
between the detailed geoid and the satellite geoid can probably be attributed to 
the following; 1) Variations of gravity are greater in Eurasia than in North 
America and are thus more difficult to detect from satellite motion, and 2) The 
surface gravity data in Asia is less accurate than that in Western Europe or 
North America. 

6. CONCLUSIONS 

The detailed gravimetric geoids presented here have a precision of ±2 
meters over the continents and 5 to 7 meters where data was sparse. 

The use of a consistent set of parameters has removed the systematic 
errors inherent in previously computed gravimetric geoids and hence reference 
these geoids to an absolute datum. 

One question that might have been answered in this study concerns the pos- 
sible existence of a rotation in different major datums. From this study there 
seems to be no conclusive evidence of a rotation in the North American datum 
but a rotation, which is prominent along the East- West profile does exist in the 
European and Australian datums. This rotation could be attributed to long wave- 
length errors in the satellite derived gravity model, a rotation of the astrogeo- 
detic geoid or a combination of both. 
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Table 1 

GSFC Long-Arc Solution*/Gravimetric Geoid Comparison (meters) 


(1) 

Station 

No. 

(2) 

Station 

Name 

(3) 

GSFC Long-Arc 
Geoid Height 9 

(4) 

Gravimetric 
Geoid Height 

(5) 

(D- © 

(6) 

(D +7m b 

(7) 

©+3m c 

United 

States 

1021 

Blossom Pt. 

-41 

-33 

-8 

-1 

-5 

1022 

Ft. Myers 

-31 

-28 

-3 


0 

1030 

Goldstone 

-35 

-33 

-2 


1 

1034 

E. Grand Fks. 

-30 

-25 

-5 


-2 

1042 

Rosman 

-41 

-30 



-8 

7036 

Edinburg 

-36 

-24 

1 



7037 

Columbia 

-42 

-32 




7045 

Denver 

-26 

-19 

mm 

0 


7050 

Greenbelt 

-40 

-32 


-1 


7072 

Jupiter 

-35 

-32 

SSI 


0 

7075 

Sudbury 

-42 

-37 



-2 

Caribbean 

7039 

Bermuda 

-41 

-40 


©- lm 
-2 

2 

7040 

San Juan 

-53 


-3 

-4 

0 

7076 

Jamaica 

-26 

-31 

5 

4 

8 

Europe 

1035 

Wi nkfield 

43 

43 

0 

(D + 0m 
0 

3 

8009 

Delft 

40 

41 

-1 

-1 

2 

8010 

Zimmerwald 

50 

47 

3 

3 

6 

8015 

Haute Provence 

54 

49 

5 

5 

8 

8019 

Nice 

47 

48 

-1 

-1 

2 

9004 

San Fernando 

47 

45 

2 

2 

5 

9091 

Dionysos 

41 

42 

-1 

-1 

2 

9115 

Oslo 

40 

37 


3 

6 

9432 

Uzhgorod 

35 

41 


-6 

-3 

Australia 

1024 

Woomera 

7 

12 

-5 

© + 6m 
1 

-2 

1038 

Orroral 

24 

25 

-1 

5 

2 

7054 

Carnarvon 

-23 

-10 

-13 

-7 

-10 

9023 

Is. Lagoon 

7 

12 

-5 

1 

-2 


a. Referenced to an ellipsoid with semimajor axis = 6378.142 km and a kM value of 3.986012 x 10 20 cm 3 /sec 2 


b. Adjusted local difference 


c. Adjusted global difference 
'Marsh, Douglas and Klosko (1971) 
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Table 2 

GSFC GEM 4 Solution*/Gravimetric Geoid Comparison (meters) 


(1) 

Station 

No. 

(2) 

Station 

Name 

(3) 

GEM 4 

Geoid Height 3 

(4) 

Gravimetric 
Geoid Height 

(5) 

(D-d) 

(6) 

(E^0 m b 


United 

States 

1021 

Blossom Pt. 

-36 

-32 

-4 

-4 

-4 

1022 

Ft. Myers 

-25 

-28 

3 

3 

3 

1030 

Goldstone 

-28 

-34 

6 

6 

6 

1034 

E. Grand Fks. 

-25 

-25 

0 

0 

0 

1042 

Rosman 

-30 

-30 

0 

o 

0 

7036 

Edinburg 

-26 

-24 

-2 

-2 

-2 

7037 

Columbia 

-33 

-32 

-1 

-1 

-1 

7045 

Denver 

-20 

-19 

-1 

-1 

-1 

7072 

Jupiter 

-35 

-32 

-3 

-3 

-3 

7075 

Sudbury 

-34 

-35 

1 

1 

1 

Carribean 

7039 

Bermuda 

-36 

-40 

4 

©+ 9m 
-5 

4 

7040 

San Juan 

-45 


5 

-4 

5 

7076 

Jamaica 

-13 

-31 

18 

9 

18 

1 

Wi nkfield 

49 

43 


(5)+ 3m 
9 

6 


San Fernando 

43 

45 


1 

-2 


Dionysos 

28 

42 

-14 

-11 

-14 

Australia 

1024 

Woomera 

12 

12 

0 

©+ 4m 
4 

0 

1038 

Orroral 

25 

25 

0 

4 

0 

7054 

Carnarvon 

-25 

-10 

-15 

-11 

-15 

9023 

Is. Lagoon 

11 

12 

-1 

3 

-1 


a. Referenced to an ellipsoid with semimajor axis = 6378.142 km and a kM value of 3.986012 x 10 20 cm 3 /sec 2 


b. Adjusted local difference 

c. Adjusted global difference 

‘Lerch, Wagner, Smith, Sandson, Brownd and Richardson, (1972) 
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Table 3 


SAO *69 Solution*/Gravimetric Geoid Comparison (meters) 


(1) 

Station 

No 

(2) 

Station 

Name 

(3) 

SAO 

Geoid Height 3 

(4) 

Gravimetric 
Geoid Height 

(5) 

(D-0 

(6) 

©+3m b 

(7) 

(D+0m c 

1 1 I 

Blossom Pt. 

-34 

-32 

-2 

1 

-2 


E. Grand Fks 

-24 

-25 



1 

1042 

Rosman 

-41 

-30 



-11 

7037 

Columbia 

-33 

-32 



-1 

7045 

Denver 

-12 

-19 




7050 

Greenbelt 

-34- 

-32 

k ms 



7075 

Sudbury 

-44 

-35 

jjE9& 



9001 

Organ Pass 

-28 

-22 




MOTH 

Jupi ter 

-32 

-32 

0 


0 

9021 

Mt. Hopkins 

-34 

-28 

-6 


-6 

9050 

Harvard 

-38 

-25 

-13 

-10 

-13 

9113 

Edwards AFB 

-30 

-35 

5 

8 

5 

Caribbean 

7039 

Bermuda 

-28 

-40 

12 


12 

7040 

San Juan 

-46 


4 

Effl 

4 

7076 

Jamaica 

-2 

-31 

29 

Wm 

29 

Austral ia 
9023 

Is. Lagoon 

12 

12 

0 


0 

9003 

Woomera 

10 

12 

-2 

D 

-1 


a. References to an ellipsoid with semimajor axis = 6378.142 km and a kM value of 3.986012 x 10 20 cm 3 /sec 2 


b. Adjusted local difference 

c. Adjusted global difference 
*Gaposchkin and Lambeck (1970) 
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Table 4 


Comparison of Detailed Gravimetric Geoid and Rice's Astrogeodetic 
Geoid Under Different Assumptions for Transforming Astrogeodetic Data 


Latitude , N 

Longitude , W 

1 

2 

3 

4 

5 

34 ° 

58 ' 

03 '.'0 

120 ° 

38 ’ 

05 '.'5 

2 

mm 

2 

3 

1 

35 

00 

38.0 

119 

00 

48.0 ■ 

2 


4 

3 

1 

38 

47 

23.1 

121 

52 

15.6 

0 



Bl /- 

-1 

35 

02 

36.1 

106 

30 

24.1 

1 

9 



1 

32 

13 

14.7 

106 

29 

41.6 

3 


99 



32 

00 

00.6 

103 

16 

07.2 

0 

0 

0 

0 


30 

59 

40.0 

098 

05 

50.5 

-2 


-1 

-2 

Warn 

30 

36 

26 ..5 

091 

23 

18.1 

-1 


-1 

-1 

mm 

29 

38 

10.8 

091 

06 

49.3 

-1 


0 

-1 

-i 

30 

59 

25.5 

089 

34 

29.5 

-2 

in- 

-2 

-2 

-2 

28 

29 

28.6 

080 

33 

35.6 

3 

6 

3 

3 

3 

30 

36 

53.3 

081 

42 

14.8 

2 

5 

2 

3 

2 

39 

28 

18.9 

076 

05 

15.2 • 

-3 

-1 

-4 

-3 

-3 

34 

59 

44.0 

076 

59 

11.7 

-3 

1 

-3 

-2 

-2 

33 

28 

42.4 

091 

00 

08.5 

-1 

1 

-1 

-1 

-1 

33 

34 

48.5 

092 

50 

07.2 

0 


1 

0 

0 

34 

56 

47.0 

093 

24 

18.3 

-2 

KoK^' 

-2 

-2 

-2 

37 

38 

08.4 

094 

35 

46.8 

-2 

If 

-1 

ISA.- 

-2 

35 

03 

04.0 

097 

56 

52.6 

0 

■fit' 

0 

■9 

0 

39 

13 

26.7 

098 

32 

30.5 

0 

99 

0 


0 

43 

37 

10.7 

096 

17 

52.3 

0 


0 


0 

35 

06 

16.2 

103 

19 

55.0 

2 


2 


1 

34 

56 

32.8 

096 

24 

55.3 

0 


0 

0 

0 

44 

43 

46.0 

105 

25 

50.7 

2 

i 

2 

1 

1 

36 

47 

44.2 

103 

11 

48.5 

1 

i 

1 

1 

0 

38 

50 

40.6 

102 

48 

46.8 

0 

0 

1 

0 

0 

48 

06 

18.6 

102 

21 

09.7 

-2 

-2 

-2 

-2 

-3 

46 

44 

47.4 

102 

15 

13.4 

-2 

-2 

-2 

-2 

-2 

45 

12 

45.7 

102 

09 

14.1 

0 

-2 

-2 

-2 

-2 

46 

21 

53.1 

108 

59 

07.3 

1 

0 

1 

1 

0 

31 

03 

07.3 

102 

56 

05.8 

0 

0 

0 

0 

-1 

41 

30 

419 . 

097 

37 

23.4 

-1 

0 

-1 

-1 

gpum 

30 

48 

49.8 

093 

12 

26.9 

-3 

-1 


-4 

Miff 

-47 

50 

28.9 

no 

00 

46.4 

1 

0 

1 

1 

9 


1 = Corrected difference between Rice's astrogeodetic geoid and the detailed gravimetric geoid using GSFC Long 
Arc (Marsh, et al., 1971) translation values of Ax = -25.1, Ay = 162.9, and Az = 172.5 meters; mean scale 
difference = 16m. 


2 = Corrected difference between Rice's astrogeodetic geoid and the detailed gravimetric geoid using Fischer's 

(1968) translation values of Ax = -18, Ay = 145, and Az = 183; mean scale difference = 2m. 

3 = Corrected difference between Rice's astrogeodetic geoid and the detailed gravimetric geoid using SAO Standard 

Earth '66 translation values of Ax = -30, Ay = 1 52, and Az = 1 76 (Lundquist and Veis, 1967); mean scale 
difference = 5m. 

4 = Corrected difference between Rice's astrogeodetic geoid and the detailed gravimetric geoid using SAO's trans- 

formation values of Ax = -25.8, Ay = 168.1, and Az = 167 (K. Lambeck, 1971, personal communication); mean 
scale difference = 23m. 

5 = Corrected difference between Rice's astrogeodetic geoid and the detailed gravimetric geoid using GSFC GEM 4 

transformation values of Ax = -24, Ay = 153, and Az = 181 (Lerch et al., 1971); mean scale difference = 3m. 
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Table 5 


Fischer’s Astrogeodetic Geoid/Gravimetric Geoid Comparison 
for Eurasia at 52 °N Latitude (meters) 














Table 6-a. 

(Detailed Gravimetric Geoid - Satellite Geoid) 
North America 


Location 

A h (meters) 

1. 47° N, 115°W 

-8 

2. 27° N, 107° W 

8 

3. 24°N, 95° W 

-10 

4. 55°N, 83° W 

-10 

5. 42° N, 72° W 

14 

6. 28° N, 71° W 

-8 

7. 17° N, 61° W 

-12 

8. 49° N, 61° W 

-10 

9. 43° N, 49° W 

12 

10. 22° N, 52° W 

-8 

11. 22° N, 30° W 

12 

12. 38° W, 20° W 

14 


Table 6-b. 

(Detailed Gravimetric Geoid - Satellite Geoid) 
Eurasia 


Location 

A h (meters) 

1. 46° N, 5°W 

-16 

2. 39° N, 3°W 

-14 

3. 41°N, 16° E 

18 

4. 41° N, 21° E 

18 

5. 54° N, 28° E 

-13 

6. 44° N, 37° E 

-18 

7. 31° N, 72°E 

-10 

8. 34° N, 79° E 

14 

9. 46° N, 86° E 

-14 

10. 26° N, 85° E 

-16 

11. 19° N, 88° E 

12 

12. 44° N, 103° E 

16 

13. 62° N, 105° E 

24 


25 












Detailed Gravimetric Geoid 

Bomford’s Astrogeodetic Geoid 
Transformed to Geocentric Using 
AX = 89, AY = 120, AZ =118 meters 


Longitude, °E 


Figure 1 - Detailed Gravimetric Geoid and Bomford’s Transformed Astrogeodetic 
Data at 44°N (Adjusted for 10m Scale Difference) 



Detailed Gravimetric Geoid 


Bomford’s Astrogeodetic Geoid 
Transformed to Geocentric Using 
AX = 89, AY = 120, AZ= 118 meters 



Figure 2 - Detailed Gravimetric Geoid and Bomford's Transformed Astrogeodetic 
Geoid at Latitude 48 °N (Adjusted for 10m Scale Difference) 
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Bomford’s Astrogeodetic Geoid 
Transformed to Geocentric Using 
AX = 89, AY = 120, and AZ= 118 
meters 



Figure 3 — Detailed Gravimetric Geoid and Bomford’s Transformed Astrogeodetic 
Geoid at 9°E (Adjusted for 10m Scale Difference) 
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Figure 4 - Detailed Gravimetric Geoid and Mather's Astrogeodetic 
Data at 22 S (Adjusted for 12 m scale difference) 
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Detailed Gravimetric Geoid 


Mather’s Astrogeodetic Geoid 
Transformed to Geocentric Using 



Figure 5 - Detailed Grgvimetric G eoid a nd Math er's Astrogeodetic 
Data at 26 S (Adjusted for 12 m scale difference) 
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Mather’s Astrogeodetic Geoid 



(Adjusted for 12 m scale difference) 
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Figure 8 - Australia Satellite Geoid 
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Figure 10 — Australia — A Detailed Gravimetric Geoid 
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